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Abstract. Pulsar Timing Arrays use a set of millisecond pulsars in an attempt to directly 
detect nanohertz gravitational waves. For this purpose, high precision timing of the pulsars is 
essential and ultimately a precision of the order of ~100ns is required. Propagation effects in 
the interstellar medium cause the radio emission from a pulsar to be dispersed and scattered, 
introducing time variable delays of the pulses on their way to Earth. If these delays are not 
properly corrected for, they may cause significant errors in the timing analysis of a pulsar. 
These proceedings will review the effects of the interstellar medium on pulse arrival times and 
present some of the techniques used to mitigate the associated time delays from the pulsar 
signal. Correcting for these delays is essential to providing a higher timing precision and hence 
to increasing the array’s sensitivity to gravitational waves. 


1. Introduction 

Pulsar Timing Arrays (PTAs) are observing a set of high-precision timing millisecond pulsars 
(MSPs) with the aim to detect nano-Hz gravitational waves (GWs). MSPs are extremely 
stably rotating neutron stars, with times of arrival (TOAs) of the pulses measurable to a few /is 
precision. In the timing procedure of a pulsar, residuals are created by making a timing model 
based on the parameters of the pulsar, and accounting for each rotation of the star. By observing 
MSPs for a long time span, the timing model is improved and the scatter of the residuals can 
go down to an RMS of ~100ns. To search for GWs, PTAs look for correlations in the timing 
residuals of several MSPs. Since PTAs are sensitive to nanohertz frequency GWs, which make 
up a part of frequency space not covered by ground based and space based GW detectors, the 
efforts of various GW projects are complementary. The main GW sources expected in the PTA 
frequency band are super massive black hole binaries and cosmic strings, and searches are being 
carried out both for single sources and for a GW background. Currently, three collaborations 
are undertaking these efforts, namely the European Pulsar Timing Array 1 (EPTA; [1]), the 
Parkes Pulsar Timing Array 2 (PPTA; [2]), and the North American Nanohertz Observatory 
for Gravitational Waves 3 (NANOGrav; [3]). Together they make up the International Pulsar 
Timing Array 4 (IPTA; [4]). 

1 http: / / www.epta.eu.org 

2 http://www.atnf.csiro.au/research/pulsar/ppta 

3 http: //nanograv.org 

http: / / www.ipta4gw.org 
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To succeed in detecting GWs, it is necessary to achieve as high timing precision as possible, 
over a long time span. As the signal from a pulsar travels through the ionized plasma of the 
interstellar medium (ISM), the pulses get scattered and the TOAs of the pulses are delayed. 
The two dominating ISM effects are dispersion and scattering, and both these effects introduce 
perturbations to the TOAs and will increase the timing error if not properly corrected for. 
Monitoring of ISM parameters and development of correction techniques is on-going within 
several pulsar groups [5-7], and ISM mitigation is increasingly recognized as one of the most 
important parts of the noise budget in PTAs [8-10]. 

This document is written as part of the 10th International LISA Symposium conference 
proceedings, and is based on a talk reviewing the challenges of ISM propagation effects on 
pulsar timing data as well as the different techniques used to mitigate these effects from the 
pulsar signals. 

2. Dispersion 

As the pulsar signal travels through the ISM, the overall column density of free electrons has 
dispersive effects on the broadband pulses. We can write the group velocity of the pulse train 
as 



where c is the speed of light, v is the observing frequency and u p is the plasma frequency given 
by 



Here, e is the electron charge, m e is the electron mass and n e is the electron density in the 
relevant direction in the sky, all given in cgs units. If looking into the Galactic plane, the value 
of n e ~ 0.03 cm -3 , and it gets lower at higher Galactic latitudes. The definition of the Dispersion 
Measure (DM) as the integrated column density of free electrons along some line of sight can be 
written as 

r d 

DM = / n e dl, (3) 

jo 

where d is the distance to the pulsar. From this, we can derive the time delay of the radio wave 
due to the DM as 

e 2 DM 

tDM( " ) = S^? DMR<4 - 15ms 4^^ (4) 

The observational consequence of the dispersive delay is that a pulsar observed with a broad 
frequency band will show pulses that arrive later at lower frequency than at higher frequency, 
following the v~ 2 law. By observing pulsars at a wide enough frequency band, or at two or more 
separate frequency bands, the DM can be measured at each observing epoch. An example of an 
observation that has not been DM corrected can be found in Fig. 1. 

One practical use of the DM is as a tool to infer the distance to the pulsar. This can be done 
with the use of models of the free electrons in the Galaxy. Since there are more free electrons 
in certain directions in the Galaxy, such as along the Galactic plane, a higher DM does not 
necessarily mean a larger distance, but also the sky position of a pulsar need to be considered. 
Currently, the NE2001 model [11] is the most comprehensive and widely used model. 

Since the pulsar and the ISM have different relative velocities, the DM of a pulsar is not 
constant in time, and the variations of the DM need to be tracked and corrected for in the data. 
This can be achieved by observing pulsars at two or more separate frequencies at a maximum 
of a few days apart. Different methods on how to optimally correct for the DM variations are 




Figure 1. Wide-band observation of PSR J1744-1134, before (left) and after (right) DM 
correction. In this case, the dispersion delay over a 400 MHz band is almost equal to one 
full rotation of the pulsar. Due to scintillation, some intensity modulation over the band is also 
visible in this observation. These plots were created from raw data, contributing to the weak 
interference signals visible in the background. 


being used by the different PTAs. Currently, the method used by NANOGrav involves fitting 
the multifrequency data for a DM within a ~10 day window, in a simultaneous fit with all 
pulsar parameters [12]. The PPTA technique is very similar to this, but removes the dispersive 
delay from the timing residuals after the timing solution has been subtracted. This method 
includes a frequency independent “common mode” delay [9], that is introduced to make sure 
that no GW signal is absorbed into the DM correction. Finally, the EPTA method uses the fact 
that the DM variations are correlated between adjacent observations, and applies a time-domain 
spectral analysis to measure the DM power-law spectrum and then extract the DM waveform 
by constructing a linear optimal filter [10]. 

3. Multi-path scattering 

Spatial inhomogeneities in the ISM give rise to scattering of the radio waves from pulsars. 
Such multi-path scattering manifests itself in several ways, including refractive and diffractive 
intensity scintillations as well as pulse broadening. 

Direct detection of pulse broadening is predominantly possible for pulsars with high DM 
values, i.e. for sources with more free electrons along their line of sight. Strong scattering cause 
large delays in the TOAs and hence high-DM pulsars are usually not included in PTAs. One 
exception is PSR J1903+0327, which because of its high flux density is useful for PTAs even 
though it has a large DM value. An observation of this pulsar can be seen in Fig. 2, showing 
a large scattering tail in its folded profile. For pulsars with more modest scattering, the pulse 
broadening can be measured indirectly through analyses of the diffractive scintillation pattern 
of the pulsar, which will be described later in this section. 

The most commonly used model to describe the ISM consists of a thin disk of turbulent 
plasma midway between the pulsar and the observer [13]. When the pulses travel through the 
disk, inhomogeneities in the plasma introduce phase perturbations to the signal. In combination 
with the different relative velocities of the pulsar, the ISM, and the Earth, this contributes 
to a diffraction pattern being observed. The perturbations are correlated over a scintillation 
bandwidth (Az/d), which is inversely proportional to the scattering delay (rd), and over a 
scintillation timescale (Aid)- Such diffractive scintillation effects were first observed in pulsars 
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Figure 2. Pulse profile of PSR J1903+0327 at an observing frequency of 1.4 GHz, showing an 
exponential scattering tail. This pulsar has a large DM = 298 pc cm -3 , contributing to strong 
scattering of the profile even at this relatively high observing frequency. 


by Lyne & Rickett [14], and originate in interstellar density structure of scale size ~ 10 6 — 10 s cm. 
The ISM plasma is often considered to be composed as a Kolmogorov medium, which implies 
scaling of the scintillation parameters as Az/ d oc z/ 4 , and Af d oc v 1 ' 2 [15]. The diffractive 
scintillation pattern for a pulsar can change drastically on short time scales of weeks to months 
[16]. Similar to the case of DM variations, the relative velocities of the pulsar and the ISM give 
rise to the time variable scattering delays. 

Larger structures in the ISM, with a size scale ~ 10 10 — 10 12 cm, give rise to refractive 
scintillation effects. These effects arise through wavefront curvature, by causing the angle of 
arrival to diverge from the straight line of sight, and are observed through flux density variations 
in pulsars [17]. 

The multi-path broadening of pulses from diffraction give rise to the largest ISM delays, next 
to dispersive delays [5]. These can be measured indirectly through dynamical spectrum analyses 
of the pulsar data. A dynamic spectrum displays the intensity of the on-pulse part of the pulse 
profile over observing time and frequency. In the spectrum, scintillation maxima show up as 
bright islands on the map (so called “scintles”, see Fig. 3). The scintles are limited in time 
since the interference pattern has a characteristic length scale, which becomes a characteristic 
time scale due to the relative velocity of the system, and they are limited in frequency because 
coherent radio beams have scattering delays of order r d . 

To get enough signal to noise to measure the scintillation parameters from the scintles, a 
2-dimensional autocorrelation function (ACF) is created from the dynamic spectrum. This 2D 
ACF is then summed in time and frequency separately, and the two resulting ID ACFs are 
fitted with a Gaussian function, centered at zero lag. From the Gaussian fit, Af d is defined as 
the half-width at e^ 1 of the summed frequency lag and Azz d is the half-width at half-maximum 
of the summed time lag [18]. The scattering delay, r d , can be calculated from At'j through 

27rAz/ d T d = C± (5) 

where C\ is a constant that depends on the geometry and spectrum of the electron density, and 
with expected value ranging from 0.6—1.5 [19]. 
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Figure 3. Dynamic spectrum of a ~7-hr long observation of PSR J1614-2230, after interference 
signals have been removed from the data (showing up as dark lines in the plot). Lighter color 
indicates higher intensity of the pulsar signal, with separate bright islands referred to as scintles. 
Analyses of the scintles in the dynamic spectrum can be used to estimate scintillation parameters 
for a pulsar. 


By creating dynamic spectra of each observation for each of the PTA pulsars, the scintillation 
parameters can be monitored, and in particular the variation of the scattering delays can be 
analyzed. This is important, since a constant delay would not cause any additional errors in 
timing precision, no matter its magnitude, but a highly variable delay would contribute largely to 
the timing noise. Therefore, the maximum variation of the delays need to be monitored, and can 
be compared to the individual TOA error to see if the scattering is one of the dominating sources 
of noise in the data. At the commonly used observing frequency band centered at ~1.4 GHz, for 
most PTA pulsars, the scattering delay variations are in general much smaller than the current 
TOA uncertainty [7, 20], and hence are currently not limiting the timing precision. With longer 
data spans and new developments in instrumentation, the timing precision will improve, and 
hence a detailed scattering analysis may become increasingly important in the future. 

The latest developments in ISM mitigation include a technique to estimate pulse broadening 
times and scattering delays using the signal-processing tool Cyclic Spectroscopy (CS; [21,22]). 
This technique allows determination of the ISM impulse response function and extraction of the 
intrinsic pulse profile for highly scattered pulsars with a large flux density. CS backends have 
been developed at some of the radio telescopes [23], providing observations with much higher 
frequency resolution than possible with the Liter banks used in conventional pulsar backends. 
Hence, very narrow scintillation bandwidths of high DM pulsars, that would not have been 
possible to measure with other pulsar backends, may be resolved in observations with a CS 
backend. This makes CS a very powerful tool for ISM mitigation and it will surely be used 
increasingly in the future. 



4. Conclusions 

PTAs strive to time MSPs with the highest possible precision, to allow for direct detection of 
nano-Hz GWs. Mitigating the delays that originate from radio wave propagation through the 
ISM is an essential part of the high-precision timing procedure. 

The largest ISM delays are due to dispersion, and time-variable DM values are currently 
corrected for in the timing data. Optimal DM correction is a highly topical issue among the 
PTA members and development of new techniques in this area is on-going. 

Further, interstellar density fluctuations also cause time-variable delays in the pulse time of 
arrival. For most PTA pulsars, scattering is currently not a limiting factor. However, with 
new wide-band receivers and higher timing accuracy, ability to correct for scattering delays 
may become an increasingly important part of high-precision timing. In addition, development 
of new techniques to properly mitigate scattering delays from the pulsar signal may allow for 
observations at lower frequencies as well as inclusion of more distant pulsars to the array in the 
future. 
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